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First principles density functional calculations within the local spin density approximation 
(LSDA) and the generalized gradient approximation (GGA) are reported for pyrochlore Cd20s207 
and Cd2Re207. The transition metal ti g manifolds are found to be well separated in energy from 
the O 2p bands and from the higher lying e 3 and Cd derived bands. The active electronic structure 
in the ti g manifold near the Fermi energy, Ep is found to be significantly modified by spin orbit 
interactions, which we include. Both materials show semi-metallic band structures, in which the Ef 
lies in an pseudogap. The band structure of Cd20s207 near Ef is dominated by very heavy hole 
and electron bands, though at Ep the electron sections are lighter. Cd2Re207 has heavy hole bands 
but moderate mass electron states. The results are discussed in terms of measured transport and 
thermodynamic properties of these compounds as well as the very different ground states of these 
two materials. 



I. INTRODUCTION 



Investigations of transition metal oxides with metal- 
insulator transitions in their phase diagrams have re- 
vealed a remarkable range of interesting, complex and 
often unanticipated phenomena, including high temper- 
ature superconductivity, various charge, spin and orbital 
ordered states, triplet superconductivity, giant magne- 
toelastic effects, and heavy Fermion behavior. A good 
example is provided by the perovskite manganites. Al- 
though known for several decades, these compounds were 
re-investigated in detail during the last few years because 
of interest in the colossal magnetoresistive effect, which 
is basically a metal-insulator transition occurring at or 
near the magnetic ordering temperature for some com- 
positions. What has been revealed is a complex phase di- 
agram including charge, spin and orbital ordered phases 
resulting from an interplay of strong correlations, strong 
lattice coupling, and band structure effects. 

The pyrochlore oxide, Cd20s207 is metallic at room 
temperature, but undergoes a metal-insulator transition 
at 226 K.M Upon cooling through the transition, the re- 
sistivity, p(T), crosses over from a metallic temperature 
dependence to an insulating (p strongly increasing as T 
is reduced) shape, though it does not fit a standard acti- 
vated form. Instead, Sleight and co-workersEl fit p(T) to 
a form consistent with a second order transition where 
the gap goes smoothly to zero at the critical temper* 
ature. Consistent with this, Mandrus and co-workersa 
report that the specific heat anomaly resembles a mean 
field transition and shows neither latent heat nor hystere- 
sis. They also reported extensive crystallographic data, 
showing very little coupling of the transition to lattice 
degrees of freedom and confirming the purely electronic 



character of the transition. The specific heat shows a 
very small but possibly still finite Sommerfeld coefficient 
(note that the reported lattice part does not fit a sim- 
ple Debye model complicating the analysis); 7 ~ 1 
mJ/mole K 2 in the range 2K to 4K and the low T sus- 
ceptibility is high for an ordered antiferromagnetic insu- 
lating phase. The thermopower above the transition is 
hole-like and small. Just below the transition, the ther- 
mopower increases to +40 fxV/K, but then decreases and 
crosses zero. A very large electron-like thermopower, jj 
to S = —300 fiV/K is seen well below the transition^ 
This is not expected from the simplest picture of a fully 
gapped Fermi surface due to a spin density wave. On 
the other hand, supposing that the metal-insulator tran- 
sition is magnetic in character, the ground state could 
be very complicated especially in view of the strong ge- 
ometric frustration of the pyrochlore lattice and the fact 
that the interactions are most likely antiferromagnetic. 
The related compound, Cd2Re2C>7, which differs from 
Cd20s207 by one electron per transition metal, was re- 
cently reported to be superconducting at ~1K with a 
substantial diamagnetic signature.Eru It displays metal- 
lic properties up to 200K where an apparently, isostruc- 
tural transition of unknown character occurs. ErEm Like 
Cd20s2C>7, this transition shows signatures in electrical 
transport and susceptibility, though the details differ. In- 
terestingly, the low temperature linear specific heat coef- 
ficient for Cd 2 Re 2 07, 7=30 mJ/molK 2 (Ref. |) is simi- 
lar to the value 7 ~ 20mJ / molK 2 estimated for the high 
temperature phase of Cd20s207 (Ref. ||) in spite of the 
difference in electron count and the very different low 
temperature properties. 



Here we report the electronic structures of these two 
compounds as calculated within the local density ap- 
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proximation (LDA) and the generalized gradient approx- 
imation (GGA), using the general potential linearized 
augmented planewave (LAPW) methpd including spin- 
orbit, which we find important .H~tL3 The calculations 
were done with well converged zone samplings and ba- 
sis sets (including local orbitals to relax linearization 
errors and treat semi-core states). The experimental 
crystal structure at 180K of Mandrus and co-workers, 
(i.e. a=10.1598A and x(Ol)=0.319) was used through- 
out for Cd20s207.B The calculated Hellman-Feynman 
force on the symmetry unconstrained structural param- 
eter was found to be small with this choice supporting 
the experimental crystallography. Total energy mini- 
mization, within the LDA, yields a value x(Ol)=0.3198. 
Donohue and co-workersa refined the crystal structure 
of Cd 2 Re 2 7 obtaining a=10.219A and x(Ol)=0.309. 
However, this yields rather short Re-0 bond lengths, 
and we note that the refinement is difficult because the 
compound contains heavy elements and the presence of 
Cd mandated the use of X-rays rather than neutrons. 
We performed a relaxation of the internal coordinate 
within the LDA keeping the lattice parameter fixed at 
its no doubt reliable experimental value. We obtained 
x(Ol)=0.316 in this way. The energy of this relaxed 
structure was 0.35 eV / cell (4 Re atoms) lower in energy 
than with the reported value of x(Ol)=0.309, which is 
well outside the normal LDA errors. The Raman active 
aig phonon frequency corresponding to this coordinate 
is calculated as 463 cm . The corresponding frequency 
for Cd20s207, also calculated in the LDA is almost the 
same, oj=459 cm -1 , indicating that the Re compound 
is not softer than the Os regardless of the unusual Re 
valence, at least as measured by this particular phonon 
mode. 



In many aspects, the calculations wc present are 
like those presented by Mandrus and co-workers for 
Cd20s20y. However, there is one important difference. 
We find that spin orbit interactions are significant be- 
cause of the presence of heavy elements in the struc- 
ture, and so include them by the usual second varia- 
tional approach. The inclusion of spin-orbit substantially 
changes the electronic structure near the Fermi energy. 
In addition, we allow for magnetism in Cd20s207, which 
we investigated using self-consistent unconstrained and 
fixed spin moment calculations. These were done by the 
method described in Ref. [L5| Not surprisingly, our cal- 
culated non-spin-polarized band structure for Cd20s207 
without including spin orbit is practically identical to that 
given by Mandrus and co-workers and so is not displayed 
here. Significantly, it is very metallic and has four Os d 
derived bands crossing the Fermi energy Ep and other 
bands. It is difficult to envisage what kind of instability 
might make it insulating. In the remainder of this report 
we discuss the band structure including spin orbit, which 
while similar in many basic aspects differs significantly 
near Ep. 



II. BAND STRUCTURE 

The crystal structure of Cd2 0s 2 07 features Os ions 
(nominally Os +5 , 5<i 3 ) at the center of O octahedra. 
Cd2Re207 is similar but with one less 5d electron per 
transition metal atom. Within an ionic model, assuming 
nominal charges, one expects a manifold of occupied O 2p 
bands, followed by a partially filled transition metal t 2g 
manifold and a higher lying unnoccupied set of e g bands. 
Furthermore, since the actual transition metal site sym- 
metry is weakly rhombohcdral (due to the second neigh- 
bor coordination) , a further crystal field splitting of the 
tig manifold is possible. 

The calculated LDA band structure of Cd20s207 is 
shown in Fig. |l|, while that of Cd2Rc207 is in Fig. ||. 
The corresponding non-spin-polarized densities of states 
(DOS) near Ep are shown in Fig. | and f§ The GGA 
band structures (not shown) are very similar. As shown 
by the band structure the t 2g and e g manifolds are cleanly 
separated from each other, and from the O 2p bands by 
clean gaps, as in the ionic model. However, despite this 
there is quite strong hybridization between Os d and O 
p states, reflecting the covalent tendency of 5d transition 
metals relative to 3d oxides. Although the actual Os site 
symmetry is rhombohcdral, this part of the crystal field is 
weak and the band structure shows no apparent further 
splitting of the t 2g derived manifold. It is exactly half- 
filled and contains 12 bands (note that there are 4 Os ions 
per unit cell). The t 2g band width is 2.85 eV. Considering 
that the effective Hubbard U is likely 2 eV or less, based 
on the trends for transition metal oxides, and noting the 
multichannel character of the t 2g manifold, Cd 2 0s 2 07 
should not be classified as a strongly correlated material 
in the sense of having on-site Coulomb correlations play 
a dominant role in the formation of the electronic struc- 
ture. This is consistent with, the observation of a high 
temperature metallic phase.EJ 

As mentioned, the t 2g manifold, which consists of 12 
bands, is exactly half-filled. Unlike the scalar relativisitic 
band structure, the band structure including spin-orbit 
shows a semi-metallic structure in the sense that there is 
a gap between the sixth and seventh bands throughout 
the Brillouin zone, but because of the dispersions it is 
not an insulating gap. Instead, there are two electron- 
like Fermi surfaces. One, from the lowest conduction 
band is a shell around the T point, while the other con- 
sists of ellipsoids along the T-X lines. Corresponding 
to these electron surfaces there are hole surfaces at the 
zone boundary around the W points. This band struc- 
ture can, at least conceptually, be made insulating in two 
ways: (1) by increasing the gap between the sixth and 
seventh bands or (2) by depressing the sixth band at the 
zone boundary and or raising the, degenerate at T, sev- 
enth and eighth bands near the zone center. The heavy 
masses of these narrow bands is consistent with high ther- 
mopowers if a gap is opened. From the band structure, 
it may be seen that the electron sheets of Fermi sur- 



2 



face are lighter (higher velocity) than the hole sections. 
For the paramagnetic metallic state, the calculated av- 
erage Fermi velocity is low reflecting these heavy bands, 
< v 2 x > x ' 2 =< v 2 y >!/2 = < v l >l/2= 5.9 x i 6 cm /s 

The calculated electronic DOS, N(E) for Cd 2 0s 2 7 
has a large peak just above Ep with N(Ep) — 9.24 eV _1 
on a per formula unit (2 Os atom) basis. The GGA yields 
a somewhat higher value N(Ep) = 11.4 eV _1 . How- 
ever, while these numbers are smaller than the scalar rel- 
ativistic value of 12.7 eV _1 reported by Mandrus and 
co-workers,cl. our LDA specific heat coefficient 7=22 
mJ/molK 2 is still consistent within the range of elec- 
tronic specific heats above the transition as estimated 
from experiment .□ This leaves little room for enhance- 
ment by electron phonon interactions or beyond density 
functional many body correlations. 

The band structure and corresponding DOS for 
Cd2Re207 are shown in Figs. || and [|, respectively. As 
may be seen from the upper panel of Fig. 0, the large en- 
ergy scale electronic structure is like that of Cd20s 2 07, 
and in particular the nominal ionic model is valid, with 
Ep falling in an isolated manifold of transition metal t 2g 
states. However, a closer examination (lower panel of Fig. 

shows that the differences from the Os compound are 
not at all well described by a rigid band model. The band 
structure, like that of the Os compound is semi-metallic. 
In particular, there is a clean pseudogap between the 
fourth and fifth bands, and since there are eight d elec- 
trons per cell the nominal Fermi energy lies between them 
The Fermi surfaces consist of nearly spherical, moderate 
mass r centered electron pockets from the fifth and sixth 
bands in the ti g manifold (m* ~ 1.2), and very heavy 
hole sections from near the zone boundary. These enclose 
a total of 0.15 e/cell and an equal number of holes. As 
may be seen, this is less symmetric than the situation in 
Cd20s20y. However, because of the semi-metallic char- 
acter, which places Ep near the band edges, the velocity 
is still rather low, < v 2 . > x / 2 = 7.5 x 10 6 cm/s. Jin and 
co-workersu report that the Hall number is quite T de- 
pendent, but is electron- like at low temperatures. This 
is consistent with our band structure, which has both 
electron and hole sheets, but with much lighter electron 
sheets, which will dominate the conductivity due to their 
higher velocities. 

The calculated N{E F ) for Cd 2 Re20 7 is 5.3 eV" 1 per 
formula unit, and is derived mainly from the heavy hole 
bands. This corresponds to a bare band specific heat 
7=12.4 mJ/molK 2 . Comparing with the measured value 
of 29.6 mJ/molK 2 (Ref. ^) one obtains an enhancement 
(1 + Xtotai) = 2.4 or Xtotai = 1.4, a reasonable value for a 
known superconductor. 

The plasma frequency for Cd2Re207 from the calcu- 
lated N(E F ) and < v 2 x > is ftw p =1.3eV. Within Boltz- 
mann tranport theory for conventional metals, the slope 
of the intrinsic resistivity, in the moderate temperature 
linear regime, is given by dp/dT — (&TT 2 hkBX tr ) / {fru p ) 2 , 
where Xt r is the transport electron phonon coupling, of- 
ten a reasonable approximation to the superconducting 



coupling A. Jin and co-workersu report approximately 
linear p vs. T from approximately 50-200K. However, 
the two reported samples differ by a factor of more than 
three in dp/dT, providing a broad range of Xt r from 0.75 
to 2.5. This is problematic, as according to the experi- 
mental data dp/dT displays unusual though fine struc- 
ture in this T range, and then changes discontinuously 
to a value near zero at the transition. It is extremely 
hard to imagine such a drastic change in X tr originating 
from the conventional mechanism, as it would suggest a 
high temperature phase with almost no electron phonon 
scattering, which crosses over into a moderate to strong 
coupled low temperature phase. 

The bare susceptibility from the value N(Ep) is 
Xb = 1-69 x 10~ 4 emu/mol. The experimental data 
show a large spread. Jin and co-workerso estimate 
x(0)=5.4x 10~ 4 emu/mol, while Hanawa and co-workersd 
obtain 3.0 xlO emu/mol at low temperature, rising to 
« 5 x 10~ 4 emu/mol above the 200K phase transition. 
This corresponds to Wilson ratios Rw from 0.74 to 1.3, 
which are low for a transition metal oxide, but would 
indicate weak electron correlations in the presence of a 
moderate electron phonon coupling. 



III. MAGNETISM IN Cd 2 0s 2 07 

It should be noted that the high value of (Np) ob- 
tained in Cd2 0s207 would lead to a strong magnetic in- 
stability in a 3c? based material, but is marginal here. 
This can be understood in terms of the expected lower 
Stoner parameter / in a 5d compound. The peak near 
Ep derives from the flat practically dispersionless band 
that lies just above Ep over most of the zone. The rel- 
atively high N(Ep) suggests the possibility of a Stoner 
instability against ferromagnetism. We checked for this 
both within the LSDA and GGA using fixed spin mo- 
ment calculations, but found that Cd20s207 is pre- 
dicted to be stable against ferromagnetism, so there is 
not such an instability. The calculated susceptibility is 
X = 9.0 x 10~ 7 emu/g in the LSDA. This was determined 
from a fourth order fit of the energy vs. moment for the 
small moment part of Fig. ^| The bare Pauli suscep- 
tibility from N(Ep) is xb = 4.2 x 10 -7 emu/g, which 
yields a Stoner enhancement of (1 — NI)^ 1 ~ 2.2. The 
GGA, which sometimes overestimates the tendency to- 
wards magnetism in Ad and 5d oxides, places Cd20s207 
closer to a magnetic instability, with x — 2.4 x 10~ 6 
emu/g and x/xb = (1 — NI)^ 1 ~ 6. Taking into account 
the difference in the LDA and GGA values of N(Ep), 
the GGA value of the effective Stoner / is 50% higher 
than the corresponding LDA value. This is reminiscent 
of the situation for S^RuO/t, where LSDA calculations 
correctly produce a paramagnetic state, while the GGA 
produces an incorrect ferromagnetic ground stateJl3'E£l 
due to an overestimated /. In any case, Stoner ferromag- 
netism would rigidly exchange split the band structure, 
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at least for small moments; as can be seen from the band 
structure, this would not result in an insulating electronic 
structure. 

The band structure is quite isotropic (note e.g. the 
similar dispersions around Ep along T-X and T-L in Fig. 
|l|) and does not display strong nesting. Thus there is no 
obvious preferred wavevector to check in searching for a 
magnetic instability. However, noting that there are four 
Os atoms in the unit cell, we checked for a T point an- 
tiferromagnetic instability, in which two of the four Os 
atoms are spin up and two are spin down. In princi- 
ple, a state like this could produce an insulating band 
structure. This is because the site dependent on-site ex- 
change splittings reduce hopping between opposite spin 
Os atoms, thus potentially narrowing the bands enough 
to open the pseudogap between the sixth and seventh 
bands producing a full band gap. While we find that 
Cd20s207 is unstable against such an antiferromagnetic 
ordering in the GGA without spin orbit, we find that 
spin-orbit favors the paramagnetic state. In the LSDA 
an antiferromagnetic ordering of this type does not occur. 
Within the GGA, including spin-orbit, we find that the 
material is on the borderline of an instability against this 
antiferromagnetic order. The resulting antiferromagnetic 
state, with spin moments of ~ 0.3 ^b/Os, has the same 
energy as the paramagnetic state to within the precision 
of our calculations. Furthermore, with this size of mo- 
ments, the band narrowing is insufficient to destroy the 
metallic state. The Os sublattice consists of corner shar- 
ing tetrahedra, with the Os atoms at the corners. As a 
result it is not possible to find a structure in which all 
nearest neighbor Os bonds are antiferromagnetic, and 
additionally, the lattice is strongly geometrically frus- 
trated for nearest neighbor spin Hamiltonians. How- 
ever, based on our band structure results, the proximity 
of Cd20s207 to magnetism is itinerant in nature, and 
therefore cannot be described in terms of simple near- 
est neighbor spin Hamiltonians. In fact, the experimen- 
tal low thermodynamic data do not show any evidence 
of a large near ground state degeneracy. However, it is 
still the case that the simple collinear antiferromagnetic 
state we considered is favored relative to a ferromagnetic 
state, and so one may speculate that more complicated 
arrangements that would better fulfill a tendency for an- 
tiferromagnetic neighbors may be lower in energy. In 
particular, even though small moment itinerant magnets 
are generally collinear, as may be understood from argu- 
ments based on band kinetic energy considerations, here 
for want of a better alternative, one might postulate that 
a non-collinear ground state may occur. One example of 
such a state would be an arrangement in which the Os 
moments are directed either towards or away from the 
center of the tetrahedra. Such a state could be at the T 
point or could have the character of a spin density wave, 
with modulation at a non-zero wavevector. 

Each Os atom has two like spin and four opposite spin 
neighbors, in the simple antiferromagnetic configuration 
we calculated. On the pyrochlore lattice, this leaves a 



fully connected 3D network of like spin nearest neighbor 
atoms. One may speculate that a more complex arrange- 
ment that disrupts this fully connected network (as the 
non-collinear arrangement mentioned above would do) 
would narrow the bands enough to produce an insulat- 
ing state. It is unclear whether such a state would be 
energetically favored here. However, the fact that we 
find the simple collinear antiferromagnetic state to be at 
least marginally stable in the GGA and favored with re- 
spect to ferromagnetism, may suggest that it is at least 
possible. In an itinerant magnetic state, like the one 
discussed above, the magnetic interactions in real space 
are inherently long range, so depression of the ordering 
temperature, due to geometric frustration is not to be 
expected. Still we emphasize that this scenario is not well 
supported by the current results for three reasons: (1) It 
is unclear that a purely magnetic instability exists in the 
LSDA, which may be a better approximation than the 
GGA for 5d compounds. An alternate scenario is a tran- 
sition involving both charge and spin ordering, though 
there is not presently evidence for a charge density wave 
in the experimental data; (2) Non-collinear states, as 
mentioned, are generally not found or expected in itin- 
erant low moment magnets; and finally (3) even if the 
system has magnetic instabilities it is unclear that they 
could be strong enough to produce a transition temper- 
ature above 200K. 



IV. SUMMARY, DISCUSSION AND 
CONCLUSIONS 

Band structure calculations for Cd20s207 and 
Cd2Re207 show a considerable sensitivity of the elec- 
tronic structure near Ep to spin orbit interactions. Both 
materials show semi-metallic band structures with heavy 
bands near Ep. Cd2Re207 has heavy hole bands near 
the zone boundary and relatively light electron pockets 
around T, while the electronic structure of Cd20s207 
is dominated by heavy bands for both the holes and 
electrons. Interestingly, because of the higher N(Ep) 
in Cd 2 0s207, the transport function, Nv 2 differs by 
only 7%. However, the specific heat enhancement 
in Cd20s207 is apparently quite small, while that in 
Cd2Re207 is 2.4, leaving little room for electron-phonon 
interactions or simple many electron effects in the former. 
Furthermore, while Cd20s207 is at least near antiferro- 
magnetism in the GGA, neither compound shows a clear 
proximity to magnetism in the LSDA, which may be a 
more reasonable approximation for bd compounds. There 
are no clear nesting features in the band structures that 
would suggest spin or charge density wave Fermi surface 
instabilities. 

It is tempting to speculate that the metal-insulator 
transition in Cd20s207 and the 200K transition in 
Cd2Re207 are related to some common feature in their 
electronic structures. In this regard, we note that al- 
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though they have different electron counts, both com- 
pounds have semimetallic band structures each with 
nominally equal numbers of holes and electrons domi- 
nated by very heavy bands. This is suggestive of a ex- 
citonic instability of the Fermi surface of the type prfe 
posed by MottEZI and reviewed by Halperin and Rice.Ea 
The theory in its simplest form involves pairing between 
electrons on one sheet of Fermi surface and holes on an- 
other. Conditions that favor such a state are (1) heavy 
band masses (2) low carrier densities and (3) similar 
sizes and shapes of Fermi surfaces, although the latter 
condition can sometimes be relaxed as discussed in Ref. 
|l8| . Such an instability is purely electronic, and in con- 
trast to e.g. a charge density wave need not be cou- 
pled significantly to lattice degrees of freedom. The ex- 
citonic state may be either singlet, or triplet.—. Further- 
more it can coexist-with superconductivityEjCj or band 
ferromagnetismJSacJ provided that there are excess car- 
riers. While such a transition is due to electron correla- 
tions, it is not associated with an on-site Hubbard-like 
Coulomb repulsion, but rather k dependent correlations 
near Ep. Importantly, the formation of an excitonic state 
involves weak thermodynamic signatures (it involves elec- 
trons only within a distance determined by the exciton 
binding energy from Ef) and is continuous. 



The general phenomenology of an excitonic metal in- 
sulator transition would seem to fit experimental knowl- 
edge and the above band stricture results for Cd20s 2 07, 
especially for triplet pairing Ej The key parameters deter- 
mining whether such a state can be formed are the exci- 
ton binding energy and the degree to which the size and 
shape of the electron and hole Fermi surfaces match in k 
as compared to the inverse size of an exciton. We cannot 
assess these as we do have the effective dielectric screen- 
ing. However, we note that the very heavy bands should 
help to produce small excitons and favor this possibil- 
ity. The situation in Cd2Re2C>7 is similar. However, in 
this material the 200K transition is between two metallic 
states. Within the excitonic scenario the metallic con- 
duction below 200K and the superconductivity depend on 
the existence of non-paired carriers below the transition. 
These could be provided by e.g. slight off-stoichiometry, 
which may perhaps be anticipated from the unusual Re 
valence. This doping could be sample dependent provid- 
ing an explanation of the very different dp/dT and the 
different T c of the two reported samples discussed above. 



We emphasize that the above discussion of excitonic 
states in Cd 2 0s207 and Cd2Re20y is speculative. How- 
ever, if it is so, these materials will be interesting novel 
tests of many body theories of excitonic phases, including 
superconductivity, possible triplet pairing and the pres- 
ence of strong spin-orbit interactions. 
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FIG. 1. LDA paramagnetic band structure of 

Cd20s207. Spin-orbit interactions were included via a 
second variational step. The horizontal reference denotes 
Ef- The bottom panel is a blow-up around Ef showing 
the Os ti a manifold. 




FIG. 2. LDA paramagnetic band structure of 

Cd2Re207. Spin-orbit interactions were included via a 
second variational step. The horizontal reference denotes 
Ef- The bottom panel is a blow-up around Ef showing 
the Re t2 g manifold. The heavy bands cross Ef yielding 
hole sections of Fermi surface, but these are not along the 
symmetry lines shown. Note the non-rigid band behavior 
relative to Fig. hi. 
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FIG. 3. LDA paramagnetic electronic density of states 
of CCI2OS2O7 for the t2 g manifold on a per unit cell (4 Os 
atoms) basis. Spin-orbit is included. 




FIG. 4. LDA paramagnetic electronic density of states 
of Cd2Re2 07 for the tig manifold on a per unit cell (4 Re 
atoms) basis. Spin-orbit is included. 
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FIG. 5. Energy (Ry/cell) as a function of magneti- 
zation (/xs/cell) for Cd2 0s2 07 from fixed spin moment 
calculations including spin orbit. The solid (dotted) lines 
are for the LSD A (GGA). 
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